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Abstract 

We present the photometric properties of a sample of infrared (IR) bright dust obscured galaxies 
(DOGs). Combining wide and deep optical images obtained with the Hyper Suprime-Cam (HSC) on 
the Subaru Telescope and all-sky mid-IR (MIR) images taken with Wide-Field Infrared Survey Explorer 
(WISE), we discovered 48 DOGs with i — Ks> 1.2 and i — [22] > 7.0, where i, Kg, and [22] represent AB 
magnitude in the Tband, Ag-band, and 22 yim, respectively, in the GAMA 14hr field (~ 9 deg^). Among 
these objects, 31 (^ 65 %) show power-law spectral energy distributions (SEDs) in the near-IR (NIR) and 
MIR regime, while the remainder show a NIR bump in their SEDs. Assuming that the redshift distribution 
for our DOGs sample is Gaussian, with mean and sigma z = 1.99 ± 0.45, we calculated their total IR 
luminosity using an empirical relation between 22 /rm luminosity and total IR luminosity. The average 
value of the total IR luminosity is (3.5 ± 1.1) x lO^^L©, which classifies them as hyper-luminous infrared 
galaxies (HyLIRGs). We also derived the total IR luminosity function (LF) and IR luminosity density 
(LD) for a flux-limited subsample of 18 DOGs with 22 flux greater than 3.0 mJy and with i-band 
magnitude brighter than 24 AB magnitude. The derived space density for this subsample is log (p = -6.59 
± 0.11 [Mpc“^]. The IR LF for DOGs including data obtained from the literature is well fitted by a 
double-power law. The derived lower limit for the IR LD for our sample is pm ^ 3.8 x 10^ [Lq Mpc“^] 
and its contributions to the total IR LD, IR LD of all ultra-luminous infrared galaxies (ULIRGs), and that 
of all DOGs are > 3 %, > 9 %, and > 15 %, respectively. 
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1. Introduction 

Recent studies have revealed that almost all massive 
galaxies harbor a supermassive black hole (SMBH) with 
a mass of at their centers. Interestingly, their 

masses are strongly correlated with those of the spheroid 
component of their host galaxies, suggesting that galaxies 
and SMBHs coevolve (e.g., Magorrian et al. 1998; Marconi 
& Hunt 2003). How did galaxies and SMBHs form and 
evolve in the 13.7-billion-year history of the universe? The 
mechanism for the co-evolution of galaxies and SMBHs 
has not been well-constrained observationally, although 
it has been the subject of intense theoretical investiga¬ 
tion (e.g., Di Matteo et al. 2005; Hopkins et al. 2008). 
This is in part because many previous studies have been 
based on optically selected samples. While deep X-ray 
(e.g.. Page et al. 2004; Alexander et al. 2005; Alexander 
et al. 2011; Ueda et al. 2014), radio (e.g., Park et al. 
2008; Smolcic et al. 2009; Bardelli et al. 2010), and mid- 
infrared (e.g.. Stern et al. 2005; Lacy et al. 2013; Lacy 
et al. 2015) studies can probe the heavily obscured galax¬ 
ies that harbor actively growing black holes which are diffi¬ 
cult to find optically, they do not cover enough sky to find 
many of the rarest, most luminous objects. For a full un¬ 
derstanding of the physics of galaxy-SMBH co-evolution, 
it is crucial to search for actively accreting galaxy-SMBH 
systems which may be surrounded by heavy dust. 

Here we focus on “Dust Obscured Galaxies” (DOGs: 
Dey et al. 2008) as a key population to tackle the mystery 
of the co-evolution. DOGs are very faint in the optical, 
but are bright in the IR. The original definition of DOGs 
was R — [24] > 7.5 where R and [24] represent AB magni¬ 
tudes in the i?-band and 24 ^m, respectively (Dey et al. 
2008; Fiore et al. 2008). Their niid-IR (MIR) flux densities 
are three orders of magnitude larger than those at optical 
wavelengths, which implies dust heating by significant star 
formation (SF), an active galactic nucleus (AGN), or both, 
and the bulk of the optical and ultraviolet (UV) emission 
from them is absorbed by dust. There are two sub-classes 
among the DOGs, characterized by their optical-IR spec¬ 
tral energy distributions (SEDs) (Dey et al. 2008). The 
so-called “power-law (PL) DOGs” show a continuous rise 
in flux density to longer wavelengths, for which IR color- 
color plots (e.g., Melbourne et al. 2012) imply the pres¬ 
ence of AGN activity in their nucleus. In contrast, the 
so-called “bump DOGs” exhibit a rest-frame 1.6 /rm flux 
excess that is probably due to the stellar photospheres of 
cooler stars, which implies that the energy source of their 
huge IR emission is dominated by SF activity. 

A possible key process for the co-evolution is a gas-rich 
galaxy merger, because the gas accumulating onto the nu¬ 
cleus triggers the quasar activity, whose energy then sig- 

* Based on data collected at the Subaru Telescope, which is 

operated by the National Astronomical Observatory of Japan 

(NAOJ). 


nificantly affects the evolution of their host galaxies (e.g., 
Hopkins et al. 2006, see also Sanders et al. 1988). In the 
context of this galaxy merger scenario, the central BHs 
and their host galaxies are obscured by a large amount of 
gas and dust during the initial stage of the co-evolution. 
Interestingly, the DOG color criterion is particularly sen¬ 
sitive to heavily reddened galaxies at z ~ 2 (e.g., Desai 
et al. 2009; Yan et al. 2013). This redshift is an inter¬ 
esting epoch because the peak of SF and AGN activity 
in the Universe and the bulk of stellar mass assembly in 
galaxies occurred around this epoch (e.g., Le Floc’h et al. 
2005; Richards et al. 2006; Goto et al. 2011; Bouwens et al. 
2011; Courteau et al. 2014; Madau & Dickinson 2014). 

A hydrodynamic simulation conducted by Narayanan 
et al. (2010) found that major mergers appear as 
PL DOGs when the accretion rate on the BH peaks. 
According to Dey et al. (2008), the fraction of PL DOGs 
increases with 24 /rm brightness; rising to almost 80% 
at 24 /rm flux > 1 mJy, although this fraction has an 
uncertainty of more than 20%. Therefore, IR-bright 
DOGs in particular could constitute a key population for 
understanding the co-evolution of galaxies and SMBHs. 
However, efficient searches for luminous DOGs have been 
difficult in previous surveys because (i) these objects are 
optically too faint to detect in optical bands, and (ii) they 
are spatially rare ((2.82 ± 0.05)x Mpc“^; Dey 

et al. 2008) because the timescale of the rapid growth of 
BHs in the growth phase is much shorter than that of the 
AGN activity. Therefore, high-sensitivity and wide-area 
surveys in both optical and MIR are required to search 
for the most IR-bright DOGs. 

In this study, we perform a systematic search for IR- 
bright DOGs based on wide and deep images in the optical 
and MIR. The optical data are obtained with the Hyper 
Suprime-Gam (HSC: Miyazaki et al. 2012), a gigantic mo¬ 
saic CCD camera with a 1.5 deg diameter FoV, which is 
mounted at the prime focus of the Subaru Telescope. An 
ambitious HSC legacy survey is on-going, which will cover 
a wide area (~ 1,400 deg^) of the sky^. The expected sen¬ 
sitivity of this multi-band survey {g, r, i, z, and j/) is r ~ 
26 in AB magnitude, almost 40 times deeper than that 
of the Sloan Digital Sky Survey (SDSS: York et al. 2000). 
The MIR data are obtained from the Wide-field Infrared 
Survey Explorer (WISE: Wright et al. 2010), which was 
launched in 2009. WISE performed an all-sky survey with 
a high sensitivity in four bands, 3.4, 4.6, 12, and 22 /rm; we 
will select DOGs using the 22 ^m data. The 5a detection 
limit at 22 /rm is better than 6 mJy (Wright et al. 2010), 

^ The HSC legacy imaging survey started in March 2014 as a 
Subaru strategic program (SSP); 300 nights have been allocated 
for 5 years in total. This HSC—SSP survey consists of three layers: 
wide, deep, and ultradeep. This study uses wide-layer data ob¬ 
tained in the 2014 spring run, as described in Section 2.1. A brief 
summary of the HSC-SSP may be found at the following web 
page: http://www.naoj.org/Projects/HSC/surveyplan.html. 
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Fig. 1. Survey area vs. 24 //m detection limit for previous 
DOGs surveys (GOODS-S: Fiore et al. 2008; GOODS-N: Pope 
et al. 2008; Penner et al. 2012; E-CDFS: Treister et al. 2009; 
CDF/E-CDF/EGS: Donley et al. 2010; COSMOS: Riguccini 
et al. 2011; Calanog et al. 2013; SWIRE (XMM/Chandra): 
Lanzuisi et al. 2009; NDFWS: Dey et al. 2008). The color 
bar represents the number of DOGs discovered by each survey. 
The star represents the survey area and detection limit of this 
work. 



Fig. 3. Survey footprint for HSC (blue) and VIKING (yel¬ 
low) data. In this study, we used HSC data in the overlap 
region with VIKING (~ 9 deg^). 

with Ho = 70 km Mpc“^, Dm= 0.3, and Da= 0.7. 


which is over an order of magnitude fainter than those 
of previous IR all-sky survey data (IRAS: Neugebauer 
et al. 1984; Beichman et al. 1988, AKARI: Murakami et al. 
2007; Ishihara et al. 2010). The combination of sensitiv¬ 
ity and large solid angle of these two surveys allow us to 
detect mamy DOGs. The survey areas and sensitivities of 
previous DOGs surveys are summarized in Figure 1. As 
shown in this figure, the full HSG survey, with WISE, oc¬ 
cupies a unique region of parameter space. In this study, 
we performed a search for IR-bright DOGs using early 
HSG data covering ^ 9 deg^, as a benchmark for future 
HSG studies of DOGs. 

It has been suggested that most DOGs represent a sub¬ 
class of high-redshift (z^ 2) optically-faint ultra-luminous 
infrared galaxies (ULIRGs) (e.g., Melbourne et al. 2012). 
Though ULIRGs are also important objects to study the 
coevolution of galaxies and SMBHs, systematic searches 
for high-z ULIRGs are difficult due to the confusion limit 
of far-IR (FIR) imaging observations. The situation has 
been improved with the advent of capable instruments 
such as the Submillimetre Common-User Bolometer Array 
2 (SCUBA2: Holland et al. 2013) but the confusion limit is 
still been severe especially for less luminous ULIRGs. The 
Atacama Large Millimeter/submillimeter Array (ALMA) 
will overcome the confusion problem even for less lumi¬ 
nous ULIRGs, although its small field of view means that 
it cannot carry out wide-angle surveys. On the other 
hand, the rarity of DOGs on the sky means that they 
are not affected by confusion, and thus they are a useful 
way to study the dusty population at z > 2 statistically. 

This paper is organized as follows. Section 2 describes 
the sample selection of IR-bright DOGs, and their basic 
properties such as their SEDs are presented in Section 3. 
In Section 4, we derive their total IR luminosity function 
(LF) and estimate their IR luminosity density (LD). The 
cosmology adopted in this paper assumes a flat universe 


2. Sample Selection 

The sample used for this study was selected from WISE 
MIR sources with optical counterparts detected by HSC. 
A flow chart of our sample selection process is shown in 
Figure 2; we found a total of 48 DOGs.^ 

2.1. HSC sample 

In this study, we utilized the HSC-SSP S14A_0 early 
data, containing the positions and photometric infor¬ 
mation of objects detected by observations in March 
and April in 2014. The HSC early data were mainly 
comprised of two catalogs; a wide layer catalog that 
includes two band photometry (i and y), and an ultra¬ 
deep layer catalog that includes five band photometry 
{g, r, i, z, and y). In particular, we used a sample of 
16,392,815 objects detected either in f-band or y-band 
in the wide layer catalog which covers the Galaxy And 
Mass Assembly (GAMA: Driver et al. 2009, 2011) 14hr 
field (G14: (RA, Dec) ~ (216°, 0?68) with a total survey 
area of ^ 18.6 deg^, see figure 3). The expected limiting 
magnitude (5 (t, 2" diameter aperture) for i- and y- band 
is approximately 26 and 24 in AB magnitude, respec¬ 
tively. The typical full width at half maximum (FWHM) 
of the point spread function (PSF) for the i- and y- band 
data in this field is ^ 0.55 and 0.68 arcsec, respectively. 
The data observed by the HSC in the S14A run were 
analyzed through an early version of the HSC pipeline 
(version 2.12.4d_hsc) using codes from the Large Synoptic 
Survey Telescope (LSST) software pipeline (Ivezic et al. 
2008; Axelrod et al. 2010), developed by the HSC 
software team. This pipeline performs CCD-by-CCD 

^ For the selection process, we employed the TOPCAT based on 

the Starlink Tables Infrastructure Library (STIL), which is an 

interactive graphical viewer and editor for tabular data (Taylor 

et al. 2005). 
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Fig. 2. Flow chart of the sample selection process. 
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reduction and calibration for astrometry and photometric 
zeropoints, mosaic-stacking which combines reduced 
CCD images into a large high S/N coadd image, and 
catalog generation for detecting and measuring sources 
on the coadd. The photometric calibration is based on 
data obtained from the Panoramic Survey Telescope and 
Rapid Response System (Pan-STARRS) 1 imaging survey 
(Schlafly et al. 2012; Tonry et al. 2012; Magnier et al. 
2013). In this study, we employed the cmodel magnitude 
to estimate i-band flux, which is a weighted combination 
of exponential and de Vaucouleurs fits to the light 
profile of each object (see Lupton et al. 2001; Abazajian 
et al. 2004). We limited ourselves to those objects 
with clean i-band photometry and removed duplicates. 
In particular, we required: (i) objects are selected in 
i-band (filterOl = “HSC-I”), (ii) they are completely 
isolated or were cleanly deblended (deblend_nchild 
= 0), (hi) they are unique objects; repeat observa¬ 
tions are removed (f lag_detect_is_tract_inner = 
“True” and flag_detect_is_patch_inner = “True”), 

(iv) none of the pixels in their footprint are in¬ 
terpolated (flag_flags_pixel_edge = “False”), 

(v) none of the central 3x3 pixels are saturated 
(flag_flags_pixel_saturated_center = “False”), (vi) 
none of the central 3x3 pixels are affected by cosmic-rays 
(flag_flags_pixel_cr_center = “False”), (vii) there are 
no bad pixels in their footprint (flag_flags_pixel_bad 
= “False”), (viii) there are no problems in measuring 
cmodel fluxes (flag_cmodeUlux_flags = “False”), 
and (ix) they have a clean measurement of the centroid 
(flag_centroid_sdss_flags = “False” ). This yields a 
sample of 6,024,570 i-selected sources over ^ 18 deg^. 
The f-band cmodel AB magnitude ranges from 16.0 to 
27.1. 

2.2. WISE sample 

WISE performed an all-sky survey at 3.4, 4.6, 12, and 
22 ^m, with the PSF FWHM of 6.1, 6.4, 6.5, and 12.0 
arcsec. In this study, we used the latest all-sky ALLWISE 
catalog (Cutri et al. 2014) that achieved significantly bet¬ 
ter sensitivity in 3.4 and 4.6 ^tm than the WISE all-sky 
data release (Wright et al. 2010) due to improved data 
processing. The 5 (t photometric sensitivity for 3.4, 4.6, 
12, and 22 /rm is better than 0.054, 0.071, 1, and 6 mJy 
(corresponding to 19.6, 19.3, 16.4, and 14.5 AB magni¬ 
tudes), respectively. Note that the number of 12 and 22 
/rm-detected sources was smaller than earlier release; the 
improved estimates of the the local background reduced 
the number of faint objects. 

We hrst created the 22 /im-selected catalog, meaning 
that we extracted >3cr-detected objects at 22 ^m. We 
then eliminated sources flagged by the saturation flag, ex¬ 
tend flag, and image artifact flag (e.g., diffraction spikes, 
scattered-light halos, or optical ghosts), leaving a sam¬ 
ple with clean photometry. In this study, we employed 
the profile-ht magnitude for each band, which ensures re¬ 
liable photometry for point-sources. The WISE catalog 
contains the Vega magnitude of each source, and we con¬ 
verted these to AB magnitude, using offset values Am 


(mAB = nivega + Am) for 3.4, 4.6, 12, and 22 ^m of 2.699, 
3.339, 5.174, and 6.620, respectively, according to the 
Explanatory Supplement to the AllWISE Data Release 
Products 

2.3. The VISTA Kilo-degree Infrared Galaxy survey 
(VIKING) sample 

The VISTA Kilo-degree Infrared Galaxy survey 
(VIKING: Arnaboldi et al. 2007) is performing a wide 
area (^ 1500 deg^) near-IR (NIR) imaging survey with five 
broadband filters {Z, V, J, H, and Ks) using the VISTA 
InfraRed Gamera (VIRCAM: Dalton et al. 2006) on the 
VISTA telescope operated by ESO. The 5cr AB photo¬ 
metric sensitivity at Z, Y, J, H, and K^ in 2" aperture 
is 23.1, 22.3, 22.1, 21.5, and 21.2, respectively. VIKING 
aims for a seeing FWHM of < 1.0 arcsec (Andrews et al. 
2014). We used here the data release 1 (DRl), consisting 
of a total of 14,773,385 sources in 226 deg^. These data 
partially overlap with our HSC data. The overlap be¬ 
tween the two covers ~ 9 deg^, as shown in Figure 3; we 
refer to it as G14-VIKING in what follows. In this study, 
we employed the default point source aperture-corrected 
magnitude (2.0 arcsec diameter) for each band. We first 
narrowed the NIR sample to the 672,843 VIKING sources 
within the G14-VIKING region. We then extracted the 
234,286 sources with signal to noise ratio (S/N) greater 
than 5 and with clean photometric flags. 

2 . 4 . Gross-identification of HSC and WISE data 

We then cross-identified the HSC sample with the WISE 
sample. However, the angular resolutions of HSC and 
WISE are significantly different, which can give rise to 
false detections. For example, a WISE source can have 
several HSC counterparts candidates lying within the 
search radius. We use the fact that DOGs have very red 
optical-NIR colors. Hence we first joined the HSC data 
with NIR data obtained from the VIKING catalog whose 
angular resolution is roughly comparable to that of the 
HSC (< 1 arcsec), and we then adopted an optical-NIR 
color cut to reject sources unlikely to be DOGs before 
cross-matching with WISE. 

First, we narrowed our HSC and WISE samples to 
sources within the G14-VIKING region, which yields a 
sample of 3,414,849 HSC sources and of 1,339 WISE 
sources, respectively. We then cross-matched the HSC 
samples with the VIKING sample. Using a matching ra¬ 
dius of 1 arcsec, 178,062 objects were cross-identified, as 
shown in Figure 4; we expect < 0.5 % of our sources are in 
correctly matched, given the VIKING surface density. For 
this “HSC-VIKING sample”, we performed a color selec¬ 
tion based on i — Ks color. Bussmann et al. (2012) inves¬ 
tigated the SEDs of DOGs with spectroscopic redshift in 
the 8.6 deg^ NDWFS Bootes field, which has the requisite 
multi-band photometry. Figure 5 shows the i — Ks distri¬ 
bution for 90 DOGs in Bussmann et al. (2012) and our 
matched HSC-VIKING sample. The color distributions 
of the two are obviously different; the colors of DOGs are 

® http://wise2.ipac.caltech.edu/docs/release/allwise/expsup/index.html 
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Fig. 4. Histograms of the angular separation of HSC sources from the VIKING (left) and WISE (right) coordinates. A search radius 
of 1 arcsec when cross-matching with VIKING and 3 arcsec for WISE was adopted. Cross-matching with the VIKING coordinates 


selected 178,062 objects within the search radius, while that 

■ , DOGs sample (Bussmann et al. 2012) ■ 

HSC-VIKING sample (xO.0005) ■ 



-2 0 2 4 6 8 

i - Ks [AB mog] 

Fig. 5. i — Ks color distribution of the sources. The blue his¬ 
togram represents our matched HSC-VIKING sample (scaled 
down by a factor of 2000). The red histogram represents 
the DOGs sample (90 objects) presented by Bussmann et al. 
(2012). The dotted line {i — Ks > 1.2) indicates the adopted 
threshold for selecting DOGs candidates. 

much redder than those of other stars and normal galax¬ 
ies. In this study, we restrict ourselves to objects with 
i — Ks< 1.2; 83,910 HSC-VIKING sources met this crite¬ 
rion. 

We then cross-matched these 83,910 objects with the 
WISE sample by using a matching radius of 3 arcsec, 
which yielded 207 objects as shown in Figure 4. The sur¬ 
face density of WISE sources in the HSC-VIKING sam¬ 
ple is approximately 224 deg“^, so adopting this search 
radius means that the probability of chance coincidence 
is less than 0.05 %. Note that 15 WISE sources have two 
HSG counterparts because the spatial resolution of the 
HSG is better than that of WISE. In this study, we chose 
the closest object as the optical counterpart. 

For these 207 HSC-WISE objects, we adopted the 
DOGs selection criterion: 

[22] >7.0, (1) 

where i and [22] represent AB magnitudes in the HSC i- 


the WISE coordinates selected 207 objects. 



i - Ks [AB mag] 


Fig. 6. Color-color diagram {i — Ks vs. i — [22]) for the 
HSC-WISE sample. Red squares represent the DOGs that 
satisfy the color selection criterion {i — [22] > 7.0). The re¬ 
maining sample is represented by yellow symbols. Numbers 
in 0 denotes the number of objects. 

band and WISE 22 pm band, respectively. This threshold 
is consistent with the original DOGs definition, R - [24] > 
7.5, as we found by fitting the PL DOGs from Melbourne 
et al. (2012) with a single power-law. Figure 6 shows the 
color-color diagram {i — Kg vs. i — [22]) for 207 HSC- 
WISE objects; 48 DOGs (hereinafter HSC-WISE DOGs) 
were selected. Figure 7 presents the 22 pm flux and 
i-band magnitude distributions for our sample. Their av¬ 
erage and median 22 pm flux densities are ~ 3.28 and ~ 
2.76 mjy, respectively, while their average and median i- 
band magnitude are ^ 22.8 and ^ 22.9, respectively. We 
visually inspected the HSC, VIKING, and WISE images 
of these 48 sources; all are clean detections with good 
photometry. 

If we had not used the prior NIR selection (i.e., i —Kg > 
1.2) for HSC objects and simply cross-identified the HSC 
with WISE data, 310 WISE sources would have more than 
one HSC counterpart. We define the “multiple fraction”, 
that is, the ratio of the number of WISE sources that 
have multiple HSC counterparts to that of WISE sources 
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Fig. 7 . (top) 22 /.im flux distribution for HSC—WISE objects 
(black line) and HSC-WISE DOGs (red line), (bottom) the 
2 -band magnitude distribution for them. The blue dotted 
line represents the 2 -band limiting magnitude for the SDSS. 
Numbers in () denotes the number of objects. 

that were cross matched within a search radius of 3 arc- 
sec. The multiple fraction without and with adopting the 
prior selection is 310/865 ^ 35.8 % and 15/207 ~ 7.2 
%, respectively, which indicates our prior selection works 
well. 

We required a VIKING counterpart in our selection, 
and thus rejected the 3,236,687 HSC sources without a 
VIKING detection. We will quantify the incompleteness 
due to this when we calculate the luminosity function 
(Section 4.2.2). 

3. Results 

3.1. Spectral Energy Distributions 

We classified the 48 DOGs (HSC-WISE DOGs) in our 
sample into two types (PL and Bump DOGs) based on 
their observed SEDs. Note that discriminating the two 
types is difficult in general because the position of the rest- 
frame 1.6 fim. bump in their SEDs depends on redshift (see 
Melbourne et al. 2012). A quantitative classification was 
originally suggested by Dey et al. (2008) who performed 
two power-law fits to the MIR flux measurements of every 
source, the first to just the Spitzer/IRAC measurements 
(3.6, 4.5, 5.8, and 8.0 ^m data) and the second to the 



Fig. 8. The distribution of Ka-band excess {Ek^) for 48 
HSC-WISE DOGs. The threshold to distinguish between PL 
and Bump DOGs is Eji^ = 3.0, which is represented by the 
dashed line. 

combined IRAC and MIPS 24 fim data. Since this SED 
discriminant is optimized for the Spitzer data, we define 
a classification method optimized for our filters (i, Z, Y, 
J, H, Ks, 3.4, 4.6, 12, and 22 fim). In this study, we 
define the Kg-band excess to allow us to discriminate PL 
and Bump DOGs: 


where is the observed Xg-band flux for each HSC- 
WISE DOG and is the ATg-band flux extrapolated 
from the power-law fit to the 4.6, 12, and 22 ^m data. 
Since our iCg-band data are good quality (S/N > 5, see 
Section 2.3), this value is a reliable indicator of bump fea¬ 
ture. Figure 8 shows the Ek^ distribution for HSC-WISE 
DOGs. In this study, we classified DOGs with Ek^ greater 
than 3 as Bump DOGs and all others as PL DOGs, where 
the threshold was derived empirically by checking each 
SED. Figure 9 shows examples of their SEDs. To check 
our ATg-band excess classification method, we performed 
two power-law fits to the MIR flux measurements of every 
source, the first to just the HSC-NIR measurements (i, 
Z, Y, J, H, ATg, 3.4, and 4.6 fim data) and the second to 
the combined HSC-NIR measurement with WISE 12 and 
22 ^m data, in the sprit of Dey et al. (2008). The type 
classification based on the ATg-band excess is consistent 
with that based on the two power-law fits for 27 objects 
(27/48 ~ 56.2 %) that show a very steep continuum or a 
clear bump (see (a) and (b) in Figure 9). On the other 
hand, our classification method is different from that of 
Dey et al. (2008) for the remaining 21 objects (21/48 ^ 
43.8 %), but our method seems to be better from a visual 
inspection (see (c) and (d) in Figure 9). Thus our method 
based on the Xg-band excess works well. 

As a result, 31 (~ 65 %) and 17 (~ 25 %) objects were 
classified as PL and Bump DOGs, respectively, which is 
roughly consistent with those of previous studies (e.g., 
Dey et al. 2008), i.e., IR-bright DOGs tend to be AGN- 
dominated. The individual and averaged SEDs for each 
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Fig. 9. Example SEDs of DOGs in our sample. Red, blue, and green lines represent the best fit linear function determined by fits 
to longer WISE measurements (4.6, 12, and 22 /rm data), HSC—NIR measurements (f, Z, y, J, ff, Ks, 3.4, and 4.6 /rm data), and 
all measurements, respectively. (a,b) Examples of the SED for two objects whose type classification based on the iCs-band excess 
are consistent with those of Dey et al. (2008). (c,d) Examples of the SED for two objects whose type classification are inconsistent 
with those of Dey et al. (2008). Note that the observed SED in (c) fits a power-law well, while that in (d) shows a bump. 
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Fig. 10. The SEDs of PL (left) and Bump (right) DOGs. The red squares represent the average of individual SEDs represented by 
yellow lines for each type. Numbers in () denotes the number of objects. 
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type are shown in Figure 10. 

3.2. DOGs color in WISE 

Figure 11 shows the WISE color-color diagram ([3.4] — 
[4.6] versus [4.6] — [12]) for HSC-WISE sources. The typ¬ 
ical MIR colors for various populations of objects are 
shown with different colors (Wright et al. 2010). About 
60 % of PL DOGs are located within the AGN wedge 
defined by Mateos et al. (2012, 2013), who suggested 
a highly complete and reliable MIR color selection cri¬ 
teria for luminous AGN candidates based on the WISE 
and wide-angle Bright Ultrahard XMM-Newton survey 
(BUXS: Mateos et al. 2012), while about 73 % of Bump 
DOGs are located outside the AGN wedge. Thus there is 
not a one-to-one correspondence between AGN classified 
using Figure 11 and using the PL/Bump criterion, and 
it is not clear which classification is more reliable. The 
WISE colors of AGNs have been well studied (e.g.. Stern 
et al. 2012; Assef et al. 2013; Toba et al. 2014). However, 
our classification is based on broad band SEDs, with up 
to 10 bands, while the WISE only uses 3 bands. 

It is interesting that the HSC-WISE objects are not 
distributed uniformly in the WISE color-color plane. The 
left edge of the distribution of the HSC-WISE objects is 
due to the i — Kg selection, thus relatively blue objects 
such as elliptical and spiral galaxies are removed by the 
i — Kg prior criterion. On the other hand, the steepness 
of the slope of the continuum in the NIR and MIR regime 
determines the right edge. This is because 22 ^m-selected 
objects with extremely steep continuum are difficult to 
detect at shorter wavelength (3.4, 4.6, and 12 /rm), and 
thus drop out of the sample. Actually, such extremely red 
objects have been discovered by WISE. They are faint 
or undetected at 3.4 /rm (Wl) and 4.6 ^m (W2) while 
remaining easily detectable at 12 and/or 22 /xm. These 
objects (so-called “W12 dropouts” or “hot DOGs”: Wu 
et al. 2012) have extremely low surface density (< 0.03 
deg“^), and our HSC-WISE sample does not yet cover 
enough sky to include such extreme red objects. 

3.3. PL DOGs fraction as a function of 22 fim flux 

As described in Section 1, some authors have reported 
that the fraction of DOGs with PL SEDs increases with 
MIR flux. If their redshifts are similar, this suggests that 
the luminous IR sources tend to be more AGN-dominated. 
This has been confirmed for ULIRGs by many authors 
(e.g.. Yuan et al. 2010; Imanishi et al. 2010; Ichikawa et al. 
2014). 

The fraction of PL DOGs as a function of 22 /xm flux 
is shown in Figure 12 and Table 1. At faint 22 /xm flux, 
we use the Melbourne et al. (2012) sample of 113 DOGs 
using deep Spitzer data in the NO AO Bootes field. We 
confirmed that the brighter 22 /xm DOGs are more AGN 
dominated (i.e., PL DOGs), although we found smaller 
fractions than did Melbourne et al. (2012). Note that the 
classification in Melbourne et al. (2012) was performed 
visually based on the SEDs, while our classification was 
conducted quantitatively. Our result based on A"-band 
excess could indicate that the PL DOGs fraction for lu- 



Fig. 12. The fraction of PL DOGs as a function of 22 
flux. The error in the fraction was estimated using bino¬ 
mial statistics (see Gehrels 1986). The results of our work 
(IR-bright DOGs) are represented by filled red squares, while 
those of DOGs previously discovered by Melbourne et al. 
(2012) are represented by open circles. 

Table 1. Fraction of PL DOGs as a function of 22 pm flux 
density. 
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This work. 

minous DOGs is not much higher than that expected from 
previous studies. 

4. Discussion 

4 . 1 . Total Infrared Luminosity 

We first created a flux-limited subsample of DOGs to 
determine the luminosity function of these sources. We 
extracted the sample with 22 /xm flux greater than 3.0 
mJy (below which the WISE 22 /xm data become < 25 % 
complete) and x-band magnitude smaller than 24.0 in AB 
mag, which yielded 18 objects. The 22 /xm flux threshold 
is responsible for most of the sample reduction; the i- 
band threshold causes only one object to be dropped (see 
Figure 7). For x-band magnitude, we also determine the 
threshold by taking the completeness into account. The 
completeness of the x-band sample is close to 100 % at the 
threshold of x ^ 24. 
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Fig. 11. WISE color-color diagram of the WISE-HSC objects. Regions with different color shading show typical MIR colors of 
different populations of objects (Wright et al. 2010). The solid lines illustrate the AGN selection wedge defined from Mateos et al. 
(2012, 2013). 


We then estimated the total IR luminosity, Lir, (8-1000 
/im) for the 18 DOGs in our flux-limited sample. We have 
no spectra of our sources, so we assume that the redshift 
distribution P{z) is Gaussian, with mean and sigma = 
1.99 ± 0.45 (Dey et al. 2008). To estimate the total IR 
luminosity, we created 1000 simulated 22 ^m luminosity 
distributions for each source based on P{z). Note that 
the MIR flux obtained from deep Spitzer data in Dey 
et al. (2008) is much fainter. So it is not a priori obvi¬ 
ous that the redshift distribution should be the same for 
our sample. We checked this using the DOGs sample in 
Melbourne et al. (2012). Melbourne et al. (2012) esti¬ 
mated Lir, based on the broad-band SEDs for 113 DOGs 
with spectroscopically measured redshifts, including the 
FIR data of Herschel Space Observatory. Their sample 
has a mean spectroscopic redshift of ^ 1.93, similar to 
Dey et al. (2008). Applying similar flux cuts (i.e., flux (24 
/im) > 3.0 mJy and i-mag < 24.0) to their sample yields 4 
objects, with an average redshift of ~ 1.79 ± 0.51, which 
is consistent with that of Dey et al. (2008) given the un¬ 
certainties. Thus the data are consistent with a redshift 
distribution of DOGs roughly independent of flux, and we 
adopt the Dey et al. (2008) distribution in what follows. 

We also used an empirical relation between observed 24 
/im and total IR luminosity presented by Melbourne et al. 
(2012). They showed that the Lir for PL DOGs is well 
predicted by 24 /im luminosity with a mean Lm/vL^ (24 
/im) = 6.5 ± 1.4, whereas the scatter for the Bump DOGs 
was much larger. We converted from 22 /im to 24 /iin flux 
assuming a power-law MIR SED; f{y) oc , where a was 


fit to the SED. We include the uncertainty in this conver¬ 
sion in the Ljr errors. The derived total infrared luminos¬ 
ity as a function of i - [22] is shown in Figure 13, where 
the average value is (3.5 ± 1.1) x lO^^L©. This result 
shows that the IR-bright DOGs discovered by HSG and 
WISE are hyper-luminous infrared galaxies (HyLIRGs). 
Hereinafter, we label these 18 DOGs as “IR luminous 
DOGs”. 

4-2. Effect of rejected objects 

4-2.1. i-band undetected DOGs 

In this study, we used an z-band selection of HSG ob¬ 
jects, and selected DOGs by adopting i - Kg > 1.2 and i 
- [22] > 7.0. However, our color cut would allow DOGs 
to be identified even if they are undetected in z-band. As 
shown in Figure 7 (b), there are no objects with i > 24.0 
although the z-band detection limit is ^ 26.0. This is 
probably because these extreme red DOGs (some of them 
could be ATg-band undetected DOGs, see Section 4.2.2) 
are very rare, the area of our survey is too small to de¬ 
tect them. Therefore, the influence of z-band undetected 
objects on our results is expected to be small. 

4-2.2. Kg-band undetected DOGs 

When we cross-matched the HSG sample with the 
VIKING sample, we rejected 3,236,687 HSG sources that 
did not lie within the 1 arcsec search radius. Gould some 
of these rejected objects (HSC-non VIKING sample) be 
DOGs? These objects would be too faint at A'g-band to 
detect with VIKING, but would satisfy the DOGs crite¬ 
rion (z — [22] > 7.0). There are 39 such objects with 22 /zm 
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Fig. 13. Total IR luminosity distribution for 18 DOGs, as¬ 
suming a redshift distribution (z = 1.99 ± 0.45) as a function 
of i - [22]. All objects can be classified as HyLIRGs. The 
vertical error bars include the redshift uncertainty and the 
uncertainty in the conversion from MIR luminosity to total 
IR luminosity. 


flux greater than 3.0 mJy and *-band magnitude brighter 
than 24.0, among the itTs-band non-detections. Therefore, 
we may be missing as much as 68 % of the DOGs due to 
our iAg-band limit. Note that the cross-matching the HSG 
and WISE catalog without using the i - Kg prior cut can 
give rise to false detections as described in Section 2.4. So, 
this incompleteness is an upper limit. Actually, 14/39 (~ 
35.9 %) of DOGs have more than one counterpart candi¬ 
date, a multiple fraction, which is consistent with that de¬ 
rived in Section 2.4. In any case, one should keep in mind 
this uncertainty when discussing the statistical properties 
of DOGs in this study. 


4-3. Total Infrared Luminosity Function 

Here, we derive the LF following the 1/I4iax method 
described by Schmidt (1968). The l/ldnax method allows 
us to compute the LF directly from the data; no parame¬ 
terization or model assumptions are needed. The volume 
density (l>{L) and its uncertainty (T 0 (l) are derived using 
the expressions: 


0(L) = 


N 

Sw 


(3) 




max ,2 


(4) 


where Idnax.i is the maximum co-moving volume that 
would be enclosed at the maximum redshift at which the 
it\i object could be detected. The sums are over all galax¬ 
ies falling into a given luminosity bin. In the context of 
the cosmology we adopt, Knax is 


Kiax(z) = TT fo —dz'dD, (5) 

^ f HoJQJzrnin ' ^QMil + z'T + QA ^ ^ 


where Da is the angular distance for a given redshift in 
our adopted cosmology, D is the solid angle of the HSC- 
VIKING region (9 deg^ ^ 0.003 str), Zmin is the lower 


Fig. 14. The total IR luminosity function for HSC-DOGs 
at z 2. Data represented by black squares are derived by 
Calanog et al. (2013). Data represented by red and blue stars 
are our data without and with taking account the 39 DOG 
candidates without iCs-band detections (see Section 4.2.2), 
respectively. The vertical error bars are calculated from the 
Poisson statistical uncertainty. The red and blue solid lines 
represent the best fit double-power law function for all data 
points without and with including the 39 DOGs without 
A's-band detections, respectively. The red and blue shaded 
regions represent the uncertainties of the resultant best fit 
function due to the large uncertainty in the IR luminosity at 
bright end. 


limit of the redshift bin considered (zmin = 1-99 - 0.45 
= 1.54 in this study), and z^ax is the maximum redshift 
at which the object could be seen given the flux limit of 
the sample, or 1.99 -I- 0.45 = 2.44, whichever is smaller. 
We calculated Zmax numerically, as described in Toba et 
al. (2013, 2014). Note that when calculating Vmax, we 
should take account the detection limit for each survey 
(HSG, VIKING, and WISE) because our sample is flux- 
limited. However, the detection limit adopted for WISE 
is much larger than those for HSG and VIKING, given 
our typical SED, we thus assume here that the volume 
for each object is determined only by the WISE detection 
limit (flux at 22 /rm = 3 mJy). 

The completeness correction function C{z) in equation 
(5) is determined following Toba et al. (2013). First, we 
formulated the dependence of WISE completeness C{f) 
for flux / obtained from the Explanatory Supplement to 
the WISE All-Sky Data Release Products^. Secondly, 
we convert flux to redshift for each object, where lumi¬ 
nosity was treated as constant for each object. While 
C{f) depends on the region of the sky (see Explanatory 
Supplement to the WISE All-Sky Data Release Products), 
it is almost uniform in the region considered in this study 
(i.e., G14-VIKING region), and we approximate it to be 
exactly uniform. 

The total IR LF (i.e., the volume density of the galaxies 
per unit luminosity range) of our HSG-WISE DOGs com¬ 
puted with the 1/Knax method, is shown in Figure 14. The 
derived space density from the IR LF for our sample (18 
DOGs) is log (j) = —6.59 ± 0.11 [Mpc“^]. Figure 14 also 
shows the LF of DOGs presented by Calanog et al. (2013) 
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http://wise2. ipac.caltech.edu/docs/release/allsky/expsup/sec6_5.html. 
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who discovered 3077 DOGs with 24 flux greater than 
0.1 mJy in the 2 deg^ of the Cosmic Evolution Survey. 
We fit the LF for all galaxies using the double-power law 
(e.g., Marshall 1987): 



The free parameters are the characteristic luminosity L*, 
the normalization factor (jj*, the bright-end slope a, and 
the faint-end slope /3, respectively. The best-fitting val¬ 
ues are summarized in Table 2, and we found that the 
shape of the LF can be well fitted by the double-power 
law form. In particular, the bright side of the LF is not 
exponential in shape, predicted e.g., by the Schechter 
(1976) function. This is probably due to the contribu¬ 
tion of AGN-dominated DOGs, reflecting the fact that 
the AGN LF can be fitted by a double-power law (e.g., 
Ueda et al. 2003; Richards et al. 2006; Assef et al. 2011). 

As mentioned in Section 4.2.2, we may be missing as 
many as 39 DOGs due to our iLg-band limit. We also 
derived the LF including all 39 DOG candidates. The 
resultant fitting values are summarized in Table 2. When 
including the 39 DOG candidates in our sample (i.e., a 
total of 57 DOGs), their volume density and bright-end 
slope are log (j) = —6.09 ± 0.06 [Mpc“^] and a = -1.9 ± 0.1, 
respectively, meaning that the volume density increases 
and the bright-end slope gets flatter. The bright end slope 
is also sensitive to the uncertainty of the IR luminosity at 
the bright end. 

4.4- Total Infrared Luminosity Density 

One of the primary purposes in computing the IR LF 
is to estimate the IR luminosity density (LD), which is 
thought to be a good tracer of SF/AGN activity hidden 
by dust. In order to estimate the IR LD contributed by 
the DOGs, we integrated the IR LF weighted by the lu¬ 
minosity, i.e., 

Pm = j LiB.(f>{L)dL. (7) 

In this study, we integrated the best-fitting double-power 
law function over the luminosity range of 10^^ < Air [Aq] 
< 10^'^. Given the uncertainties in the IR LF at the bright 
end and the completeness of our sample, we simply esti¬ 
mated a lower limit to the IR LD of (IR luminous) DOGs. 
Figure 15 shows the resultant IR LD (pir) for DOGs. 
The value and its contribution to other populations (see 
below) are summarized in Table 3. For comparison, we 
also show the total IR LD (i.e., before subtracting out 
an estimate of the AGN contribution), piR for AGNs and 
ULIRGs derived from Murphy et al. (2011). The derived 
lower limit of IR LD for DOGs, piR (DOGs), is ~ 2.5 x 
10 ® [AqMpc-®], in good agreement with that estimated 
from Dey et al. (2008). We find that DOGs contribute at 
least ~ 20 % and ^ 61 % of the total IR LD and ULIRGs 
IR LD, respectively, which is roughly consistent with the 
comparison between the DOGs sample of Dey et al. (2008) 
and the 24 pm-selected sample of Gaputi et al. (2007). 


Note that we simply conducted a linear interpolation for 
data (Murphy et al. 2011) when estimating the total and 
ULIRGs’ IR LD at z = 1.99. 

In this study, we also calculated the IR LF for IR lu¬ 
minous DOGs (Air > 10^® Aq, i.e., the objects in our 
sample), and thus estimated their contribution to the IR 
LD of all DOGs and ULIRGs. The resultant IR lumi¬ 
nous DOGs LD of lower limit is ~ 3.8 x 10^ [Lq Mpc“®] 
and its contribution to piR (total), piR (ULIRGs), and 
piR (DOGs) are > 3 %, > 9 %, and > 15 %, respectively. 
The derived IR LD contributes less to the total IR LD at 
z ^ 2 while their contribution to IRLD of DOGs is still 
moderately large. 

5. Summary 

Using early HSC-wide survey data and the WISE MIR 
all-sky survey data, we performed a search for IR-bright 
DOGs in the GAMA I4hr-VIKING region (~ 9 deg^). 
We first created clean subsamples of HSG and WISE 
data, and also created a NIR selected subsample obtained 
from VIKING. We cross-identified HSG with the NIR sub¬ 
sample using a threshold inferred from those of previ¬ 
ously discovered DOGs to avoid mis-identification when 
cross-matching HSG with WISE directly. We then cross- 
identified with WISE data and adopted the DOGs color 
selection (i — [22] > 7.0), which yielded 48 DOGs. For 
those DOGs, we investigated their photometric properties 
as well as their statistical properties by constructing the 
IR LF and by estimating their IR LD. The main results 
are as follows: 

1. Among 48 DOGs, 31 DOGs (- 65 %) are PL DOGs 
and 17 DOGs (^ 25 %) are bump DOGs, according 
to their SEDs (Section 3.1). The WISE colors for 
PL DOGs indicate that they harbor AGNs in their 
nucleus (Section 3.2). 

2. The fraction of PL DOGs increases with increasing 
22 pm luminosity (Section 3.3). 

3. Assuming that the redshift distribution for our 
DOGs sample is Gaussian, with mean and sigma 
z = 1.99 ± 0.45, the average total IR luminosity of 
them is (3.5 ± 1.1) x 10 ^®Aq. Thus these objects 
are recognized as HyLIRGs (Section 4.1). 

4. The derived space density from the IR LF for our 
sample (18 DOGs) is log 4> = —6.59 ± 0.11 [Mpc“®], 
and the IR LF for DOGs including data obtained 
from the literature is well fitted by a double-power 
law with Aju = (3.1 ± 0.8) xlO^^ AQand bright-end 
slope a = -2.6 ± 0.2 (Section 4.3). 

5. The derived lower limit of IR LD for our sample is ^ 
3.8 X 10^ [AqMpc”®], and its contributions to the 
total IR LD, IR LD of all ULIRGs, and that of all 
DOGs are > 3 %, > 9 %, and > 15 %, respectively 
(Section 4.4). 

The current sample is based on only 9 deg^ of imaging 
data. The HSG survey will cover more than 100 times 
as much sky, 1400 deg^, when it is complete, allowing 
the identification of roughly 3300 DOGs. This sample 
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Table 2. Best double-power law fit parameters for the total IR LF of DOGs. 


(j) [Mpc“'^dex“^J 

[-^ 0 ] 

a (bright-end) 

P (faint-end) 

18 IR luminous DDGs (1.3 ± 0.6) x 10“'^ 

(3.1 ±0.8) X 10^^ 

-2.6±0.2 

-0.1±0.2 

57 IR luminous DOGs ’ (3.1 ± 1.9) x 10“'^ 

(1.2 ±0.5) X 10^2 

-1.9±0.1 

0.5±0.7 


Including the 39 possible DOGs that could be missed in the sample selection (see Section 4.2.2). 
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Fig. 15. The lower limit of infrared luminosity density (IR LD) for HSG-DOGs. Blue points represent total IR LD obtained from 
Murphy et al. (2011). Purple and red shaded regions represent the IR LD contributed by AGNs and ULIRGs, respectively. The IR 
LD for DOGs of all luminosities including our sample is shown as a red square, and the IR LF for IR luminous DOGs is shown as a 
red star. The green square shows the IR LD obtained from Dey et al. (2008) 


Table 3. The lower limit of the contribution of the IR LD for DOGs and IR luminous DOGs to other population. 



PIR 

[Lq Mpc-3] 

PiR (DOGs) 

contribution [%] 
piR (ULIRGs) 

piR (total) 

DOGs 

~ 2.5 X 10“ 

— 

- 61 

~ 20 

IR luminous DOGs 

~ 3.8 X 10' 

- 15 

~ 9 

~ 3 
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will be large enough to carry out a variety of statisti¬ 
cal analyses to understand the physical nature of these 
objects. In particular, cross-correlation with the BOSS 
spectroscopic quasar sample (Paris et al. 2014) will allow 
a proper determination of the redshift distribution of the 
sample, and the five-color optical photometry will give in¬ 
dependent constraints on the redshift. In fact, our objects 
are bright enough to allow spectroscopy to be performed. 
We will also search for DOGs in the deep fields (27 deg^) 
of the HSC survey: going a magnitude deeper in optical 
filters, it will allow identification of even redder sources 
than are explored in this paper. 
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